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Abstract 

Type IV shock-wave-interference heating on a 
blunt body causes extremely intense heating over a 
very localized region of the body. This paper presents 
an analytical solution to a heat-transfer problem that 
approximates the shock-wave- interference heating of 
an engine cowl leading edge of the National Aero- 
Space Plane. The problem uses a simplified geom- 
etry to represent the leading edge. An analytical 
solution is developed that provides a means for ap- 
proximating maximum temperature differences be- 
tween the outer- and inner-surface temperatures of 
the leading edge. The solution is computationally 
efficient and, as a result, is well suited for concep- 
tual and preliminary design or trade studies. Tran- 
sient and steady-state analyses are conducted, and 
results obtained from the analytical solution are 
compared with results of two-dimensional thermal- 
finite-element analyses over a wide range of design 
parameters. Isotropic materials as well as laminated 
composite materials are studied. Results of para- 
metric studies are presented to indicate the effects of 
the thickness of the cowl leading edge and the width 
of the region affected by the shock-wave-interference 
heating on the thermal response of the leading edge. 
Finally, a nondimensional temperature parameter is 
developed that is useful in evaluating the effects of 
several design parameters on the thermal response of 
the leading edge. 

Introduction 

One of the critical design problems for reusable 
air-breathing hypersonic vehicles, such as the Na- 
tional Aero-Space Plane (NASP), is caused by the se- 
vere shock-wave-intcrference heating that can occur 
on the leading edge of the engine cowl. This heating 
can occur when the oblique vehicle-nose bow shock 
wave or the oblique engine-inlet-ramp shock waves in- 
tersect the cowl-leading-edge bow shock wave (ref. 1). 
This shock- wave intersection is likely to occur be- 
cause the oblique shock waves, which compress the 
airflow external to the engine, must be located close 
to the cowl leading edge to maximize mass capture 
at the engine inlet to achieve efficient engine perfor- 
mance. The effects of shock- wave-interference heat- 
ing must be considered in the design of cowl-leading- 
edge concepts for hypersonic vehicles such as NASP 
because the intense local heating results in very high 
local temperatures, very high temperature differences 
through the cowl-leading-edge thickness, and very 
high thermal stresses in the leading edge (refs. 2 to 4). 

Recent experimental results for dynamic pres- 
sures lower than those expected in hypersonic 


flight (refs. 5 to 7) have been compared with 
computational-fluid-dynamics (CFD) predictions 
(refs. 2 and 8 to 10). The CFD predictions include 
the calculation of structural temperatures and ther- 
mal stresses with a finite-element method to solve a 
nonlinear, coupled system of aerothermal-structural 
equations (ref. 2). Results of this analysis substan- 
tiate the potential severity of the thermal stresses 
in the cowl leading edge. However, a CFD analysis 
is computationally expensive and not well suited for 
conceptual and preliminary design studies. Simpler, 
faster solution procedures are needed for the prelimi- 
nary design of thermal-structural concepts for the en- 
gine cowl leading edges of future hypersonic vehicles. 

This paper presents an analytical solution to a 
transient, linear heat-transfer problem that simulates 
two-dimensional (2-D) heat conduction in the skin 
of a .11 engine cowl leading edge subjected to intense 
local surface heating representative of Type IV 
shock- wave-interference heating. The heat-transfer 
problem uses a simplified geometry to represent the 
leading edge, and the analytical solution to this prob- 
lem provides a means for approximating the max- 
imum temperatures in the skin of the cowl lead- 
ing edge. Results obtained from this solution are 
compared with results of a finite-element analysis of 
the problem in which a more realistic cowl-leading- 
edge geometry is modeled. Analytical and finite- 
element solutions to a transient heat-transfer prob- 
lem in which the region of intense heating sweeps 
across the cowl leading edge are used to determine 
the effect of the sweep speed of the heated region 
on the maximum temperature difference through the 
cowl-leading-edge thickness. Results obtained from 
analytical and finite-element solutions to a steady- 
state problem in which the region of intense heat- 
ing remains stationary are also compared. Results of 
parametric studies are presented for various materi- 
als, including composites with graphite fibers that 
have a very high thermal conductivity along the 
length of the fibers. The parameters varied are the 
cowl-leading-edge thickness and the width of the re- 
gion of intense heating. Finally, a nondimensional 
temperature parameter is developed that aids in the 
investigation of the effects of several design parame- 
ters on the maximum temperature difference in the 
cowl leading edge. 

Symbols 

A cross-sectional area, in 2 

A n , B n coefficients used in analytical 

solution (see eqs. (A32) and (A33)) 

a n , b n Fourier coefficients 



r square root of ratio of orthotropic 

thermal conductivities (see eq. (3)) 

Cp specific heat at constant pressure, 

Btu/lb-°F 

D, H analytical solution parameters, in -1 

h cowl- leading-edge thickness, in. 

i unit vector 

k thermal conductivity, Btu/ft-hr-°F 

L length used in analytical solution 

(see eq. (3)), in. 

M Mach number 

q magnitude of external aerodynamic 

heating on cowl leading edge, 
Btu/ft 2 -sec 

q () magnitude of square heat pulse, 

Btu/ft 2 -sec 

q n , f undisturbed stagnation point heat 

flux, Btu/ft 2 -sec 

r \ , r '2 inner and outer cowl- leading-edge 

radius, respectively, in. 

T temperature, °F 

t. time, sec 

V square-heat-pulsc sweep speed, 

in /sec 

w square-heat-pulse half-width, in. 

X , Y separation functions for analytical 

solution (see eq. (A7)) 

x. y fixed Cartesian coordinate system 

for analytical solution 

li angular location from centerline of 

cowl, deg (see fig. 3) 

AT temperature difference, °F 

0 temperature difference, 

m,r})-T xn T F 

X n separation constants 

v percent-fiber-volume content in 

composite materials 

£. T) moving Cartesian coordinate system 

for analytical solution 

p density, lb/ ft * 

nondimensional temperature 
parameter 


Subscripts: 


A 

analytical 

FE 

finite element 

in 

inner surface 

FT 

longitudinal and transverse direc- 
tion, respectively 

max 

maximum 

2 

outer surface 


Problem Definition 

One of the critical issues for reusable air-breathing 
hypersonic vehicles, such as NASP, is the severity 
of the aerodynamic heating of the engine cowl lead- 
ing edge. In addition to the undisturbed aerody- 
namic heating of the cowl leading edge that results 
from hypersonic flight, shock-wave-interference heat- 
ing causes very high local increases in aerodynamic 
heating. This severe shock- wave- interference heating 
occurs when the oblique vehicle-nose bow shock wave 
or the oblique engine-inlet-ramp shock waves inter- 
sect the cowl- leading-edge bow shock wave, as shown 
in figure 1. 

When this shock- wave intersection occurs, six 
types of interference patterns can develop that af- 
fect external aerodynamic heating at various inter- 
action locations on the cowl (ref. 1). Of these types 
of patterns, the Type IV interference pattern pro- 
duces the most severe heating because it results in 
a supersonic jet that impinges on the cowl leading 
edge. A schematic diagram of the Type IV inter- 
ference pattern is shown in figure 2. The details of 
the flow in the supersonic jet are described in refer- 
ence 6. If the supersonic jet impinges nearly normal 
to the cowl- leading-edge surface, it causes a local- 
ized jet stagnation region on the leading-edge surface, 
and the resulting aerodynamic heating may be more 
than 30 times greater than the undisturbed stagna- 
tion value (ref. 7). 

A possible scenario for the formation of a Type IV 
shock-wave-interference pattern on an engine cowl 
leading edge of a hypersonic vehicle is described in 
reference 2 and is illustrated in figure 3. The hyper- 
sonic vehicle is assumed to be accelerating at Mach 16 
with an oblique vehicle-nose bow shock wave initially 
located outboard of a scramjet engine (fig. 3(a)). In 
figure 3(a), (i is a measure of the angle (in degrees) 
from the centerline of the leading edge and is positive 
in the counterclockwise direction. As the vehicle ac- 
celerates through Mach 16, the oblique vehicle-nose 
bow shock wave sweeps across the cowl leading edge 
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and intersects the cowl- leading-edge bow shock wave; 
thus, local aerodynamic heating is increased from the 
effects of shock- wave- interference heating. At a cer- 
tain intersecting location, the Type IV interference 
pattern is formed and a supersonic jet impinges on 
the leading edge (fig. 3(b)). The jet then dissipates 
as the vehicle- nose bow shock wave moves farther 
inboard (fig. 3(c)) and a less severe Type III shear- 
layer-interference pattern is formed. 

The change in the local heating distribution 
caused by the different shock-wave-interference pat- 
terns is illustrated in figure 3(d). In the figure, q 
is the magnitude of the external aerodynamic heat- 
ing on the cowl leading edge and t is the time at 
which the magnitude of the heating is measured. 
The dashed lines are instantaneous heat rate dis- 
tributions, and the solid bell-shaped curve indicates 
the envelope of the heat rate distribution that is as- 
sumed in reference 2. The maximum value q () at 
approximately ft = -22° (at t, = 0.14 sec) cor- 
responds to the Type IV shock- wave- interference 
heating. Such severe local heating can cause very 
high outer-surface temperatures, very high tem- 
perature differences through the cowl- leading-edge 
thickness, and very high thermal stresses near the 
shock-wave impingement point on the leading edge. 
This shock-wave-interfercnce heating complicates the 
thermal and structural analysis of the cowl lead- 
ing edge because an accurate calculation of temper- 
atures through the cowl-leading-edge thickness can 
require an integrated fluid-thermal-structural analy- 
sis (ref. 2). An integrated analysis may be required 
because of the coupling between the external fluid 
and thermal analyses. A simpler method for obtain- 
ing an approximate solution of temperatures within 
the skin of the cowl leading edge is desirable for use 
in preliminary design studies. 

In the present paper, the shock- wave- interference 
heating problem is assumed to be two-dimensional 
in a plane perpendicular to the spanwise direction 
(fig. 1). Other simplifying assumptions include a sim- 
plified geometry to represent the cowl leading edge, 
a square- heat- pulse representation of the supersonic- 
jet heating on the outer cowl surface, orthotropic 
temperature-independent material properties, a con- 
stant temperature on the internal cowl surface (T[ n ), 
and negligible surface radiation heat transfer. The 
simplified problem is solved through use of classic 
techniques as well as a finite-element analysis that, 
models a more realistic cowl-leading-edge geometry. 

In a previous experimental study of shock-wave- 
interference heating on a cylindrical leading edge, 
the width of the region heated by the shock-wave- 
interference heating (referred to herein as the jet 


width) was estimated to be one-twentieth of the di- 
ameter of the cylinder used in the experiment (ref. 3). 
However, this estimate was based on experimental 
results at Mach 6 and 8 and may not be represen- 
tative of values for higher Mach numbers. Errors in 
predicting the jet width can cause large errors in cal- 
culations of the thermal response of the leading edge. 
The effect, of changes in the jet width on the max- 
imum temperature difference in the leading edge is 
therefore determined subsequently by a parametric 
study. 

An experimentally determined heat-rate distribu- 
tion on a cylinder subjected to Type 1 IV shock- wave- 
interference heating is given in reference 2 and is 
shown in figure 4(a). In this figure, the ratio q/q n >{ 
is equal to the magnitude of the Type IV shock- 
wave- interference heating normalized to the magni- 
tude of the undisturbed stagnation point heat flux 
and ft is the angular location on the leading edge. 
For the present study, this distribution is approxi- 
mated by a uniform square heat pulse of constant 
width 2 w. For the transient problem (fig. 4(b)), the 
square heat pulse moves across the cowl loading edge 
at a constant speed V . For the steady-state problem 
(fig. 4(c)), the square heat pulse is stationary and 
centered at. ft — 0°. Since shock- wave-interference 
heating dominates the thermal response in t he cowl 
leading edge, it is studied independently from the' 
undisturbed stagnation heating. 

The baseline cowl- leading-edge geometry selected 
for study is shown in figure 5. The cowl consists 
of a wedge with a cylindrical nose having an inner 
radius of 0.1 in. and straight sections 2 in. long in- 
clined at angles of d=8° to the horizontal. The unit 
vectors iy’ and iy shown in the figure represent the 
directions of the orthotropic thermal conductivities. 
The thermal conductivity in the transverse direction 
(radial direction as shown in fig. 5) is A’y • and the 
thermal conductivity in the longitudinal direction is 
kj . The inner-surface temperature T\ n is assumed 
to remain constant. The assumption of a constant 
inner-surface temperature is made as an approxima- 
tion for a liquid-metal heat-pipe-cooled concept or for 
a convcctively cooled concept using a coolant with a 
high heat-transfer coefficient. However, it is, in gen- 
eral, a lioneonservative assumption for convcctively 
cooled concepts. 


Materials Investigated 

Both metallic and composite materials are con- 
sidered in this study. The metallic materials investi- 
gated are copper and Incoloy Alloy 909 (a superallov 
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made by Inco Alloys International, Inc.). The prop- 
erties used for these materials are given in table 1(a). 
Copper and Incoloy 909 are used to represent upper 
and lower bounds, respectively, of thermal conduc- 
tivity for useful candidate materials for NASP. It is 
assumed that results from the steady-state analysis 
for materials with thermal conductivities within this 
range will be bracketed by these upper and lower 
bounds. 

The composite materials investigated in this 
study consist of laminates of copper and Incoloy 909 
with 50-percent fiber-volume fraction of high- 
conductivity graphite, such as highly oriented py- 
rolytic graphite (HOPG), or carbon fiber, such as 
Thornel P-140 carbon fiber. The HOPG (proper- 
ties shown in table 1(a)) is not available as a fiber, 
but its material properties are used in this study 
as an upper bound to illustrate the maximum ca- 
pability that might be achievable for high-thermal- 
conductivity graphite fibers (ref. 11): The Thornel 
P-140 carbon fiber (properties shown in table 1(a)) is 
a commercially available, ultrahigh-modulus, pitch- 
based carbon fiber with high longitudinal thermal 
conductivity (about twice that of copper) and rela- 
tively low transverse conductivity (ref. 12). Material 
properties for the composites are calculated based 
on parallel and series conduction paths through the 
material and fiber (ref. 13). Longitudinal and trans- 
verse thermal conductivities are calculated from the 
following relationships: 


k L = uk FL + (1 - v) k M (1) 


, > _ kM^FT 

1 [( 1-0 kpj' + vkfrf] 


( 2 ) 


In equations (1) and (2), ki and kx are the longi- 
tudinal and transverse thermal conductivities of the 
composite material, respectively, kpp and kpx arc 
the longitudinal and transverse conductivities of the 
fiber, respectively, k\,f is the conductivity of the ma- 
trix, and v is the percent- fiber- volume content of the 
composite material. Bulk composite densities and 
specific heats are calculated by averaging the values 
of the constituent materials. The material properties 
for the composite materials used in this study are 
given in table 1(b). These properties are assumed to 
remain constant with temperature. 


2-D Finite-Element Analysis 

The finite-element analysis system used in this 
study is the Engineering Analysis Language (EAL, 
ref. 14). The thermal-analyzer portion of EAL 


(ref. 15) is used to perform the thermal- finite-element 
analyses. For nonlinear problems (i.e., those in- 
cluding surface radiation), the analyzer iterates on 
nodal temperatures using a modified Newton method 
until convergence of the temperatures is obtained. 
Thermal-finite-element models of both the transient 
and the steady-state problem were developed, and 
both of these models, as well as the assumed bound- 
ary conditions, arc described subsequently. 


Transient Finite-Element Analysis 

The finite-element analysis of the transient prob- 
lem is performed with an implicit, transient-response 
thermal analyzer entitled TRTB. This analyzer in- 
tegrates in time using the Crank-Nicolson method, 
and it uses the modified Newton method to solve for 
nodal temperatures at each time step (ref. 15). 

A detailed view of the cylindrical portion of 
the cowl-leading-edge thermal-finite-element model 
is shown in figure 6. The inner-surface tempera- 
ture T[ n I s assumed to remain constant, and material 
properties are assumed to remain constant with tem- 
perature. Heat is rejected from the external surface 
by radiation to space; a surface emissivity of 0.8 is as- 
sumed. Adiabatic boundary conditions are assumed 
at the ends of the model, as indicated in figure 6. The 
wall thickness h is equal to the outer radius minus 
the inner radius ( 7*2 — n), and the temperature dif- 
ference of interest in this study AT ma x is equal to the 
maximum outer-surface temperature minus T lu . The 
supersonic jet is approximated by a square heat pulse 
of width 2w. This square heat pulse is assumed to 
sweep across the entire curved portion of the leading 
edge at a constant speed V for the transient problem. 
The initial condition is T(t = 0) = T\ n . The model 
of the cowl leading edge has 2688 degrees of free- 
dom and consists of 2544 quadrilateral conduction 
elements, 384 line heat source elements, and 424 line 
radiation elements. 

The discretization of the model shown in figure 6 
is very fine around the cylindrical portion of the 
cowl leading edge. This discretization is used for 
two reasons: (1) to accurately represent a moving 
supersonic jet for the transient cases and (2) to model 
accurately 2-D heat conduction from the narrow 
region heated by the square heat pulse. The outer 
arc length of each element in the cylindrical portion 
of the leading edge is selected to be 0.2 w. At time 
t = 0, the square heat pulse is located over the first 
10 elements of the fine mesh on the lower surface 
of the model. Movement of the square heat pulse 
around the leading edge to simulate the surface effect 
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of the supersonic-jet heating is achieved by shifting 
the square heat pulse clockwise one element at each 
time step. The time step used is calculated by 
dividing the arc length of the outer surface of the 
leading edge by the number of elements in the leading 
edge and the sweep speed of the supersonic jet. The 
square heat pulse is assumed to have a constant 
magnitude, width, and speed as it traverses the entire 
curved portion of the leading edge. 

As illustrated in figure 3(d), the assumption of a 
uniform square-heat-pulse speed and a uniform mag- 
nitude of the heating is conservative compared with 
the actual heating environment. The actual transient 
heating reaches a peak when a Type IV supersonic-jet 
impingement occurs normal to the surface and then 
decreases as the incident shock continues to move 
across the leading edge. In addition, the uniform 
speed used in the present study will only equal the 
actual supersonic-jet speed at the symmetry point of 
the cowl. At all other locations, the actual speed 
would be higher. 


Steady-State Finite-Element Analysis 

The finite-element analysis of the steady-state 
problem is performed with a steady-state thermal an- 
alyzer entitled SSTA. Nonlinear problems (e.g., those 
including surface radiation) are solved as previously 
discussed. For the steady-state analysis, only half of 
the model shown in figure 6 is used to take advantage 
of the physical and thermal symmetries. Symmetry 
conditions at a location of (3 — 0° on the leading edge 
are imposed by an adiabatic wall condition. Addi- 
tionally, for the steady-state case, the square heat 
pulse is assumed to remain stationary at a location 
of (3 = 0°. The width of this pulse is equal to w , 
and it is referred to herein as the square- heat- pulse 
half-width. 

Analytical Solution of Simplified 
Problem 

Transient Solution 

The curved geometry of the cowl leading edge 
shown in figure 6 is approximated by the planar rect- 
angular region shown in figure 7. Neglecting curva- 
ture when calculating heat conduction in a curved 
cylinder is considered adequate when the ratio of the 
cylinder wall thickness h to the inner radius of the 
cylinder r\ is less than 0.1. For this study, r\ re- 
mains constant at 0.1 in. Thus, wall thicknesses less 
than or equal to 0.01 in. are considered small, and a 
rectangular approximation of a curved leading edge 


is considered appropriate. For the purposes of evalu- 
ation, wall thicknesses up to 0.03 in. are used in the 
present paper. It is expected that neglecting curva- 
ture of the cowl leading edge will be less appropriate 
as the value of h increases and approaches 0.03 in. 
The supersonic- jet heating is represented by a local 
square heat pulse that moves along the upper sur- 
face of the rectangular region at a constant speed V . 
Therefore, the simplified model of the shock-wave- 
interference heating problem consists of the rectan- 
gular region heated on the surface y = h by a moving 
square heat pulse of magnitude q a and width 2 w. 

A coordinate system (£, r/) moves with the square 
heat pulse as described in the appendix. The rect- 
angular region is subject to constant-temperature 
boundary conditions along the surfaces at £ = L, 
£ = — F, and 7/ = 0. These boundary conditions are 
applied to simplify the solution of the problem: how- 
ever, erroneous results may occur if \L\ is too small. 
Therefore, L is retained as a parameter in the so- 
lution so that appropriate conditions for it can be 
determined. Radiation heat transfer from the high- 
temperature gas of the supersonic jet to the cowl- 
leading-edge surface is neglected because, for a max- 
imum cowl-leading-edgc outer radius of 0.13 in., the 
radiation-heat-transfer contribution can be assumed 
to be a small portion of the total cowl-leading-edge 
heat load (ref. 16). Radiation leaving the upper sur- 
face of the rectangular region is also neglected. 

An analytical solution to this problem is ob- 
tained that provides an approximation to the tem- 
peratures in the cowl leading edge as a function of 
several parameters: square-heat-pulse sweep speed 
V, cow 1- leading-edge wall thickness h , square-heat- 
pulse half-width w, longitudinal thermal conductiv- 
ity hi, transverse thermal conductivity k T , and ex- 
ternal heating rate q Q . 

The development of the transient analytical so- 
lution to the simplified problem is presented in the 
appendix. In the analysis of this problem, the ef- 
fects of any initial temperature distribution are as- 
sumed to be negligible (refs. 17 and 18). The solu- 
tion to this problem as presented in the appendix is 
a quasi-steady-state solution since it further assumes 
that the temperature distribution is constant within 
the reference frame attached to the moving square 
heat pulse. However, in the present paper, this solu- 
tion is referred to as the transient analytical solution, 
since it represents a solution to the moving-square- 
heat-pulse problem. As shown in the appendix, the 
transient analytical solution in the moving coordi- 
nate system (£, rj) is 
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the finite-element and analytical solutions shows that 
the analytical solution is accurate for a wide range of 
parametric variables and material selections. This 
range includes potential NASP cowl-leading-edge ge- 
ometries and materials. 

Steady-State Solution 

To obtain an analytical solution to the steady- 
state problem in which the square heat pulse is 
stationary, the square- heat- pulse sweep speed V is 
set equal to zero in the transient analytical solution 
(eq. (3)). The analytical solution then reduces to 
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By evaluating this equation at £ = 0 and 7] = h , we 
may obtain a prediction of the maximum difference 
between the outer- and inner-surface temperatures. 
The 1 solution evaluated at this point is 


0({),fi) ATjiulx 
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(4) 


The boundary conditions assumed at £ = ±L, 
the rectangular approximation of a curved region, 
and the neglect of surface radiation heat transfer 
may all affect the accuracy of the solution. How- 
ever, as shown subsequently, a comparison between 


By evaluating this equation at £ = 0 and 7] = /i, we 
may obtain a prediction of the maximum difference 
between the outer- and inner-surface temperatures 
for the steady-state case. The solution evaluated at 
this point is 


0 (0, h) ~ AThiax 
_ SLqo 

T&C.k'l' 
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x E 
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‘sin (n7rw/2L) tanli (nnch/2L) 
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(6) 


Results and Discussion 

In the present paper, the analytical solutions 
to the shock-wavc-interference heating problem are 
compared with the finite-element solutions for both 
the transient and steady-state problems. The com- 
parisons show the applicability of the analytical 
solutions over a wide range of the following param- 
eters: square-heat-pulse sweep speed, cowl-leading- 
edge thickness, square-heat-pulse width, material 
orthotropic thermal conductivities, and magnitude 
of the shock- wave-interference heating. Parametric 
studies using the steady-state analytical solution are 
then conducted to investigate more closely the effects 
of changes in the cowl-leading-edge thickness and 
square-heat-pulse width on the maximum tempera- 
ture differences. Finally, a nondimensional temper- 
ature parameter 3> is derived from the steady-state 
solution. This parameter is useful in evaluating the 
effects of all the aforementioned parameters on the 
maximum temperature difference. 
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Transient Finite-Element and Analytical 

Results 

The effect of changes in the square- heat- pulse 
sweep speed on the maximum differences between 
the outer- and inner-surface temperatures is investi- 
gated with the transient finite-element analysis and 
the transient analytical solution (eq. (4)). Maxi- 
mum temperature differences AT max are calculated 
for a 0.01-in. -wide square heat pulse with a magni- 
tude of 50 000 Btu/ft 2 -sec sweeping across a 0. 02-in. - 
thick cowl leading edge at uniform speeds ranging 
from 0 to 20.0 in/sec. (See table II.) The values of 
A T max through the cowl- leading-edge thickness arc 
calculated for a constant inner-surface temperature 
of 0°F, and they are shown as a function of square- 
heat-pulse sweep speed in figures 8 and 9. Results for 
copper, P- 140/copper, and HOPG/copper are shown 
in figure 8, and results for Incolov 909, P-140/Incoloy 
909, and HOPG/Incoloy 909 are shown in figure 9. 
In these figures and those that follow, the maximum 
temperatures shown may exceed the allowable tem- 
peratures for the materials indicated in some circum- 
stances. However, for the purposes of discussion and 
evaluating the analytical solutions, these results are 
retained unless otherwise noted. 

As shown in figures 8 and 9, the results from 
the analytical solution overpredict AT max (compared 
with the results from the finite-element analysis) for 
all the materials and square-heat-pulse sweep speeds 
except for copper and Incoloy 909. The analytical 
and finite-element results agree for the copper, and 
for Incoloy 909 the analytical solution ovcrpredicts 
AT max at V = 0 in/sec and underpredicts AT max 
at V > 2 in/sec. However, the analytical and finite- 
element solutions give the same trends and agree (for 
the worst case of Incoloy 909) to within 15 percent. 

The range of square- heat- pulse sweep speeds 
shown in figures 8 and 9 spans the range of expected 
NASP shock- wave sweep speeds (0.2 in/sec being the 
current expected value). Within the speed range 
of 0 to 2 in/sec (10 times the expected square-heat- 
pulse sweep speed), only Incoloy 909 shows an appre- 
ciable decrease in AT max for increasing speed. This 
decrease in AT max occurs because a smaller square- 
heat-pulse residence time reduces the amount of heat 
input locally into the leading edge, and this reduc- 
tion compensates for the poor ability of Incoloy 909 
to conduct heat from the square-heat-pulse region. 
Within this speed range, the maximum cowl-leading- 
edge temperatures for all the other materials reach 
a constant value almost instantaneously, and this 
value is very close to the steady-state (zero-square- 
heat-pulse speed) predictions of AT max shown in fig- 
ures 8 and 9. Therefore, the remaining studies use 


the steady-state analytical solution to investigate all 
the materials except Incoloy 909. 

Steady-State Finite-Element and 

Analytical Results 

The steady-state analytical and finite-element so- 
lutions are compared to assess the accuracy of the 
analytical solution. The analytical solution is used 
in parametric studies to investigate the effects of 
changes in the cowl- leading-edge thickness h (with 
the inner radius held constant) and the square-heat- 
pulse half- width w on the values of AT max . 

Accuracy of the analytical solution. Pre- 
dictions of AT max obtained from the steady-state 
analytical solution and from the steady-state finite- 
element analysis are compared in figure 10 to as- 
sess the accuracy of the analytical solution. The 
data represent results obtained from the two solu- 
tions while the following parameters are varied: lon- 
gitudinal and transverse thermal conductivities kj 
and cowl-leading-edge thickness h with the inner 
radius held constant, square-heat-pulse half- width u\ 
heating rate q () , and internal operating temperature 
T\ n . The variations in thermal conductivity are ob- 
tained by changing materials, and the ranges of pa- 
rameter variations used are shown in table II. 

Results from the finite-element analysis are plot- 
ted on the ordinate in figure 10, and results from t he 
analytical solution are plotted on the abscissa. The 
dashed line represents a perfect correlation between 
the two solutions, and the actual correlation between 
the two solutions, represented by the symbols, devi- 
ates very little from this line. The maximum differ- 
ence (with respect to the analytical results) between 
the values of AT max obtained from the analytical 
and finite-element solutions is less than 10 percent. 
Therefore, based on the assumption that the finite- 
element solution gives an accurate representation of 
the temperature distribution, such a correlation be- 
tween the analytical and finite-element results indi- 
cates that the analytical solution is effective in pre- 
dicting the maximum temperature differences in the 
cowl leading edge for a wide range of parameters. 

Effect of changes in the cowl-leading- edge 
thickness on maximum temperature differ- 
ences. The steady-state analytical solution is used 
to study the effect of changes in the cowl- leading-edge 
thickness h on the values of AT niax . Results from 
this aspect of the study are presented in figure 1 1 for 
copper, P-140/copper, and HOPG/copper and in fig- 
ure 12 for P-140/Incoloy 909 and HOPG/Incoloy 909. 
The values of AT max through the cowl leading edge 
are divided by the magnitude of the heat flux q () so 
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that any value of heat flux may be considered. The 
values of A T m8LX /q 0 are shown in figures 11 and 12 
as a function of the cowl- leading-edge thickness h. 
For the results shown in these figures, the square- 
heat-pulse half-width w is assumed to be constant at 
w = 0.005 in. 

As expected, the values of AT max /q 0 increase as 
the cowl-leading-edge thickness increases because the 
length of the transverse conduction path increases. 
Furthermore, the slopes of the curves are steeper for 
leading-edge thicknesses from 0.0025 to 0.0075 in. 
compared with the slopes for larger thicknesses. This 
effect occurs because the heat conduction in the lead- 
ing edge is nearly one-dimensional in the transverse 
direction when the leading-edge thickness is small. 
Therefore, the maximum temperature difference is 
more sensitive to changes in the thickness. As the 
leading-edge thickness increases, the slopes of the 
curves decrease because the heat conduction becomes 
two-dimensional . 

The results shown in figures 11 and 12 may also 
be used to determine the maximum temperature dif- 
ferences for potential cowl concepts. For example, if 
the maximum allowable cowl-leading-edge thickness 
h is 0.005 in. for a copper leading edge subjected to 
q D = 50 000 Btu/ft 2 -sec, the value of A T max /q a is 
equal to l°F-in 2 -sec/Btu. This value corresponds to 
a maximum temperature difference of approximately 
347°F. (See fig. 11.) 

The effects of two-dimensional heat conduction in 
the cowl leading edge are also illustrated in figures 11 
and 12. The values of A T m8iX /q 0 for P-140/copper 
and HOPG/copper are larger than those for copper, 
as shown in figure 11. This result occurs because 
the transverse thermal conductivities of the compos- 
ites are low compared with the thermal conductiv- 
ity of copper. (See table I.) Furthermore, the values 
of A T maLX /q a for the HOPG/copper are lower than 
those for P- 140/copper because the HOPG/copper 
has a higher longitudinal thermal conductivity than 
the P-140/copper and has the same transverse ther- 
mal conductivity. Similarly, the values of A T m&x /q 0 
for HOPG/Incoloy 909 are lower than those for 
P-140/Incoloy 909 because the longitudinal thermal 
conductivity of HOPG/Incoloy 909 is higher than 
that of P-140/Incoloy 909. 

Effect of changes in the square-heat-pulse 
width on maximum temperature differences . 
The width of the surface region affected by the 
shock-wave-interference supersonic jet used in the 
present study (i.e., the square- heat- pulse width) is 
taken from experimental results at M = 6 and 8 
(ref. 3), and it may not be representative of values 


for higher Mach numbers. Therefore, the effect of 
changes in the square-heat-pulse width on maximum 
temperature differences predicted by the analytical 
solution is studied. 

Values of A T max /q 0 for copper and P-140/copper 
are shown in figure 13 as a function of the square- 
heat-pulse half-width w for three cowl- leading-edge 
thicknesses. Similar results for P-140/Incoloy 909 
and HOPG/Incoloy 909 are shown in figure 14. The 
sensitivity of the values of A T max /q 0 to changes 
in the square-heat-pulse half-width is equal to the 
slope d(AT max /q 0 )/dw of the curves in these fig- 
ures. The values of the slopes are given in the fig- 
ures for the curves corresponding to a cowl-leading- 
edge thickness of 0.03 in., with the assumption 
that w = 0.005 in. As shown in figure 13, the slope 
of the curve for a copper leading edge is approx- 
imately equal to 260°F-in 2 -sec/Btu-in. Thus, the 
value of A T max /q a changes 260°F-in 2 -sec/Btu per 
1-in. change in w. The effect of this sensitivity to 
changes in the square-heat-pulse half-width is il- 
lustrated in the following example. For a copper 
leading edge, if it is assumed that w — 0.005 in., 
h = 0.03 in., and q Q = 50 000 Btu/ft 2 -sec, a value of 
AT max /< 7 0 = 1 .67°F-in 2 -sec/Btu (ATi nax = 580°F) is 
calculated with the steady-state analytical solution. 
However, if the actual value of w is 50 percent 
higher (i.e., w = 0.0075 in.), the value of A T max /q 0 
increases to 2.32°F-in 2 -sec/Btu (AT max — 805°F). 
This new value represents a 38-percent increase 
over the original value. Similarly, the sensitivity 
of the P-1 40/copper composite is approximately 
1590°F-in 2 -sec/Btu-in. Therefore, for the same con- 
ditions described previously, the value of A T max /q 0 
would increase by 40 percent. These high sensitivities 
of the maximum temperature difference to changes in 
the square-heat-pulse width indicate that the width 
of the square heat pulse can have an important effect 
on the thermal response of the leading edge. 

The effects of two-dimensional heat conduction 
discussed previously are illustrated in figures 13 and 
14 as well. A large decrease in transverse thermal 
conductivity (as in the case of P-140/copper 
compared with copper) or longitudinal thermal con- 
ductivity (as in the case of P-140/Incoloy 909 com- 
pared with HOPG/Incoloy 909) increases the sen- 
sitivity of the maximum temperature differences to 
changes in square-heat-pulse width. Another trend 
illustrated in figures 13 and 14 is that the sensitivity 
of the maximum temperature difference to changes 
in the square-heat-pulse width increases as the cowl- 
leading-edge thickness increases. This increase in the 
sensitivity occurs because the length of the transverse 
conduction path increases. 
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Nondimensional Temperature Parameter 

The steady-state analytical solution (eq. (6)) is 
rearranged to obtain a nondimensional temperature 
parameter that, for L/w = 80 (see appendix), is 
given as 


, AT^iax 

<J> = 

wq 0 /ck T 

640 ^ fsin (n7r/160) tanh (n7rr/i/160iy) 

7T 2 [ n 2 

n= 1,3,... 

The only term that will affect the summation on the 
right-hand side of equation (7) is the ratio hc/w. 
Thus, only wq 0 /ck^ and hc/w govern the values of 
AT max . The dependence of on hc/w is shown in 
figure 15. The values of $ presented in the figure 
are the converged values of the infinite series on 
the right-hand side of equation (7). (See appendix 
for a discussion of convergence.) The maximum 
temperature difference A T max can be determined 
from figure 15 for a given set of parameters. For 
example, if hc/w = 1,<£ is approximately equal to 
3.47. The value for AT max can then be determined 
by multiplying 3.47 by the value of wq () /ckj. 



Concluding Remarks 

An analytical solution to a heat-transfer prob- 
lem that approximates the problem of shock- wave- 
interference heating of the engine cowl leading edge 
of the National Aero-Space Plane (NASP) has been 
presented. The problem uses a simplified geometry 
to represent the leading edge and represents the su- 
personic jet that is characteristic of Type IV shock- 
wave-interference heating as a square heat pulse. The 
analytical solution provides a means for approximat- 
ing maximum temperature differences between the 
outer and inner surfaces of the leading edge. Re- 
sults of transient and steady-state analyses from the 
analytical solution were compared with results of 
two-dimensional finite-element analyses over numer- 
ous combinations of the following design parame- 
ters: square-heat-pulse sweep speed V (which varied 
from 0 to 20.0 in/sec), cowl-leading-edge thickness h 
(which varied from 0.005 to 0.030 in.), square- heat- 
pulse width 2w (which varied from 0.01 to 0.03 in.), 
transverse (radial) thermal conductivity kj (which 
varied from 7 to 226 Btu/ft-hr-°F), longitudinal (nor- 
mal to transverse) thermal conductivity ki (which 
varied from 9 to 690 Btu/ft-hr-°F), and heat flux q 0 
(which varied from 10 000 to 100 000 Btu/ft 2 -sec). 


Results obtained with the transient, finite- 
element, and analytical solutions for square-heat- 
pulse sweep speeds expected for the NASP vehicle in- 
dicated that most of the materials used in this study 
responded nearly instantaneously to the shock-wave- 
interference heating. Thus, steady-state predictions 
of maximum temperature difference A T max were very 
close to the transient results. 

An assessment of the accuracy of the steady-state 
analytical solution indicated that the analytical so- 
lution predicted maximum temperature differences 
to within 10 percent of the finite-element results. 
Parametric studies using the steady-state analytical 
solution were also performed to investigate the ef- 
fects of changes in the cowl-leading-edge thickness 
and square- heat- pulse width on the maximum tem- 
perature differences for two different groups of ma- 
terials. The materials investigated included copper, 
Incoloy 909, and laminated composites of these mate- 
rials consisting of 50-percent fiber-volume fraction of 
either highly oriented pyrolytic graphite (HOPG) or 
Thornel P-140 carbon fibers. Copper and Incoloy 909 
were chosen for study to represent upper and lower 
bounds, respectively, for thermal conductivity of use- 
ful NASP vehicle candidate materials for the engine 
cowl leading edge. 

Results of the parametric studies illustrated that 
maximum temperature differences through the thick- 
ness of the leading edge increased when the cowl- 
leading-edge thickness was increased. Furthermore, 
the rate of change of the maximum temperature dif- 
ference with increasing leading-edge thickness was 
larger when 0.0025 < h < 0.0075 in. than it was when 
h > 0.0075 in. These studies also illustrated the ef- 
fects of two-dimensional heat conduction in the lead- 
ing edge by showing that large decreases in either 
the transverse thermal conductivity or the longitudi- 
nal conductivity of the leading-edge material resulted 
in increased maximum temperature differences. 


The value of A7j nax was very sensitive to changes 
in the width of the square heat pulse used to simulate 
the shock- wave- interference heating. In one case, a 
50-percent increase in this width resulted in an ap- 
proximately 40-percent increase in AT max . Addition- 
ally, large decreases in either the transverse or the 
longitudinal thermal conductivity of the leading-edge 
material resulted in increased sensitivity. Thus, the 
sensitivity of maximum temperature differences to 
changes in the square-heat-pulse width was affected 
by the two-dimensional nature of the heat conduction 
in the leading edge. The sensitivity of maximum tem- 
perature differences to changes in the width of the 
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square heat pulse also increased as the cowl-leading- 
edge thickness increased. 

Finally, a nondimensional temperature parameter 
was derived from the steady-state analytical solution. 
Use of this parameter makes it possible to quickly 
calculate maximum temperatures and temperature 


differences in the cowl leading edge for the range of 
parameters considered in the present study. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
January 21, 1992 
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Appendix 


Development of the Transient Analytical Solution 

The transient. Type IV shock-wave-interfcrcnce heating of the engine cowl leading edge shown 
in figure 6 is approximated by the transient, linear, two-dimensional heat-conduction problem 
represented in figure 7. The geometry is approximated by a planar rectangular region heated on the 
surface y — h by a square heat pulse of magnitude q Q and width 2 w that, moves at a constant speed 
V. This appendix presents the development of the analytical solution to this problem. 


Derivation of Governing Equation 

The governing heat-transfer equation for transient, two-dimensional heat conduction is given b\ 


dT , _ 

= *'• a? 


q2 t () 2 T 

+ k T 


dy 1 


(Al) 


where k L and k T are the material thermal conductivities in the x- and ^-directions, respectively. A 
reference frame may be attached (as in refs. 17 and 18) to the moving square heat pulse by replacing 


x, y , and t with new coordinates £, r/, and r such that 

£ = x -Vt T] = y t = t 


If T = T(£, T], t ), 

dT dTde, dTdji dZ^Z = ^y Sl Z + ^K 

~dt ~ d£ dt + dr] dt + dr dt d£ dr 

and equation (Al) may be rewritten as 


dT 

PC *Tt 


OT , d 2 T d 2 T 

VpCp dS ~ kl W + kT dv 2 


(A2) 


If it is assumed that the temperature transient due to an initial condition has sufficiently decayed to 
have a negligible effect at the time of interest, the temperature distribution may be considered to be 
steady in the moving reference frame. (See refs. 17 and 18.) Therefore, the derivative with respect 
to t in equation (A2) is zero, and the following equation results: 


dT , d 2 T d 2 T 


(A3) 


or 


c 2 % + F7> + T m = 0 


(A4) 


where 

2 _ 

C *T 


r _ v P c p 
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Equation (A4) is the governing heat-transfer equation of the transient, linear heat-conduction 
problem shown in figure 7. This equation can be solved subject to the following boundary conditions: 

= Tin 

T (L, T), t) = Tj n 
T (£> 0, r) = Tj n 


j Vo {—w < £ < w) 

’l= h l 0 (141 > w) 

Solution Procedure 

Equation (A4) is a partial differential equation that may be solved in several ways. One solution 
procedure is as follows. 

A new temperature variable may first be introduced through the following substitution: 

0(40/) = r(£,r?) -T iu (A5) 

where T in is the internal-surface temperature (at r/ = 0) that, remains constant. Substituting 
equation (A5) into equation (A4) yields 

t>2 % + Fe i + = o (A6) 

The corresponding boundary conditions are 


7(4) 


OT 

On 


9 (-L, rj , r) = 0 
8(L,t),t) = 0 
9 (4, 0, r) = 0 


q (4) = kj< 


09 

dr} 


?f—h 


j <io 
~ 1 0 


( — w < £ < w) 

(141 > w) 


Now, through use of the method of separation of variables (ref. 19), a solution is assumed in the 
form 


0(4,7) = *(4) Y(r}) (A7) 

where A - (4) and Y (rj) are separation functions. Substituting this expression for 9 into equation (A6) 
yields the following two ordinary differential equations: 

c 2 X^ + FX^ + A 2 c 2 X = 0 (A8) 

~ A 2 c 2 Y = 0 (A9) 

where A is a separation constant. 

Equation (A8) may be solved by assuming the solution: 


X(£) = exp (D£) 


(A10) 
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Substituting this assumed solution into equation (A8) yields the following equation: 

D 2 + 2 HD + A 2 - 0 


o u — — f Jl — i* 

In — 2 — 2 “ 7 , 

vr k r f( kj 

Solving for D yields 

D= -H± \/H 2 - A 2 (A 12) 

Therefore, the solution to equation (A8) is 

X(£) = A exp | + V // H 2 - A 2 ) 4- B exp [(~^ - \/ H 2 ~ ^ 2 ) (A13) 

where A and B are constants. The boundary conditions #(£ — — L) = 0(£ = L) = 0 may be used to 
determine A. Applying these boundary conditions results in 

A exp [(-# + V H 2 - A 2 ) (-L)J + B exp ^-// - \J H 2 - A 2 ^ ( — L)J = 0 (AM) 

A exp [(-// + VH 2 - A 2 ) L] + B exp [(-// - s/W 2 - A 2 ) l] = 0 (A15) 

To obtain a nontrivial solution, the determinant of the coefficients for A and B in equations (A 14) 
and (A15) is set equal to zero, the result being a linear dependence of the two equations. The 
following equation, which is solved for A, results: 


exp ^—2 LyJ H 2 - A 2 ^ — exp ^2 L\f H 2 — A 2 j 


If |A| < \H\, then the only solution for A is A — H: therefore, no solution exists. If |A > |//|. thru 
the exponent is imaginary and, with Euler's identity, equation (A 16) may be expressed as 


2i sin 2L v A 2 - H 2 = 0 


Solving for A yields 

Rewriting the expression for D yields 


,2 _ n1 * 2 , „2 
A " - JJT + 11 


D = -H ± v 7 // 2 - A 2 = -H ± —i 

£ L 

This expression may be substituted into either equation (A 14) or equation (A 15) since these two 
equations are linearly dependent. Substituting the expression for D into equation (A 14) yields 

A n exp[(-i7 + i £,] + B n t;xp[(-// - i ^) ij = 0 (A18) 

Only the ratio of B n /A n can be determined from equation (A 18). Therefore, for A n = 1, it can be 
shown that B n is 

B n = — cxp(m7r) — — cos nn — i sin titt 

For n = 2, 4, sin nn = 0 and cosn7r = 1. Therefore, B n = — 1 — ?’(0) = —1, and it may be shown 
that 


X n (0 = exp(-B£) sin ^ 


[n = 2, 4, ...) 
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For n — 1,3, ...,sin rnr = 0 and cosn 7 r = -1. Therefore, B n = 1 - z(0) = 1, and it may be shown 
that 


An(0 =exp(-J/£)cos^ (n = 1, 3, ...) (A20) 

Thus, equations (A 19) and (A 20 ) represent the solution for equation (A 8 ). 

Equation (A9) may be solved with the assumption of a solution of the form 

Yn {*}) = E n sinh \ n crj + F n cosh \ n cq (A 21 ) 

where E n and F n are constant coefficients. With the boundary condition F(0) = 0, it can be shown 
that Y n ( 7 /), the solution to equation (A9), is 

Y n ( 7 ]) = E n sinh A n c7] (n = 1,2,...) (A22) 

Finally, with the substitution of equations (A19), (A 20 ), and (A 22 ) into equation (A 7 ), the 
temperature difference between the temperature at any point in the moving coordinate system and 
Tin is 


0 fori) = T for,) -Tfr 


= exp(-H£) 


C oc 

^ A n sin — - sinh X n crj 

n=2,4,... ZL 


DC 

+ E Bn 

n=l,3,... 


cos T -^y- sinh A n cr] 

ZLj 


(A23) 


The coefficients A n and B n may be evaluated with the nonhomogeneous boundary condition 


where, at 77 = h, 


«(0 = *r 


80 

dr\ 


_ I Qo 
l=h \ 0 


(—in < £ < in) 

(l£l > w ) 


d0 / it jj £ 

— = c exp (-//£) V A n A n sin -—f- cosh A n ch 

dr > L-w,... 2i 


B n X n cos — y- cosh A n ch 


n=l,3,... 


(A24) 


Therefore 


/ n7r £ 

q (0 exp(H£) = fc r c V] A n X n sin — y- cosh X n ch 
\n=2,4,... ZL 


00 

+ B n X n COS cosh A n ch 


n=l,3,... 


(A25) 
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Equation (A25) is an infinite Fourier series. To solve for A n and B n , the Fourier coefficients a n and 
b n must be determined. 


The Fourier coefficients are defined as follows: 

a ” = i/l /K)cos ^‘ i? 


(n = 1,3,...) 


(A26) 




where 


f{0-9 (0 exp(tf£) = exp (HO 


(n = 2,4,...) 
( — w < £ < w) 

(l£l > w ) 


(A27) 


a n = Bn/c^AnCCOsh A n c/i 

6 n = Anfe^AnCCOsh A n C^l 


(A28) 

(A29) 


Therefore, by carrying out the integration indicated in equations (A26) and (A27), a n and b n may 
be shown to be 

n f f nnw . Tr , —nnw 

“» - r + ( {" r p( "“’ )cos it - cxp( - H " )c “ 


n-7T ,,, x . mrw , TT \ • ~ nnw 

+ — exp(Hw) sin — exp (-Hw) sin 

2 L 2L ^ 


(n — 1, 3, ...) (A30) 


L [H 2 4- (n 2 7r 2 /4L 2 )] 


j// exp(Hw) si 


n7Tte / rr \ — h- 7 ™ 

sin — exp(-//w) sin — 


nn , „ , n7rie , rr ^ 

— expf/ftc) cos — exp(-//w)cos — — 

2L 2L ^ 


= 2,4,...) (A31) 


By substitution of equations (A30) and (A31) into equations (A28) and (A29) and simplification of 
the results, it can be shown that A n and B n are 


Lkj'cXn cosh X n ch [ H 2 + (rc 2 7r 2 /4L 2 )] 


nn / nnw . 

— 2 cos si 
2 L V 2L 


sinh Hw'j J 


/ n7nc . __ \ 

Hi 2 sin - cosh Hw J 


(n = 2,4, ...) 


B = S ^ — \h(2 cos ^ sinh Hu;) 

n y LkjcXn cosh X n ch [H 2 + (n 2 7r 2 /4L 2 )] L V 2L ' 

+ ^ (2 sin - 7 ^ cosh Hw^ \ ( n = 3, •*•) 
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Finally, substituting equations (A32) and (A33) into equation (A23) gives the temperature distri- 
bution in the moving coordinate system (£,77): 


*(£*/) = 


q () exp (-//£) 
Ljk'j'd 

nn / 

- — (2 cos 


l= fj { A„[fl2+(!,V/4i2)] [" (2 «») 


2L 


/i7rw; 

~2L 


sinh Hw^J ( sin 


nn£ sinh A n C7/ ' 
2L cosh A n ch / 


+ E < 

n=l ,3,... 


A„ [ff 2 + (nV/4L 2 )] 


K 


n7rre \ 

2 cos — — sinh i/ic J 


+ 


n7r 

2L 


^2 si 


7/70/7 

sin - ■ cosh //u; 


)] 


cos 


nn£ sinh \ n cr) ' 
2 L cosh A nch . 


(A34) 


where 

6(Z,ri) = T(t,T))-T hl = AT(S,r,) 

In the present study, the maximum temperature difference AT max is assumed to be located at £ = 0 
and 7] — h ; however, it will actually shift to £ < 0 as the speed of the square' heat pulse increases. 
The transient solution evaluated at £ = 0 and t] = h is 


0 (0, h) = AT max 


<lo j 

f oo 

E 

Kn =1,3, 

Lk T c \ 

n 7i 
+ 2 L 

^2 sin 


A„ [// 2 + ( 


1 r / mrw \ 

n^JuJ)] [ H 1 2C0S ~2L Sinh HW ) 


nnw 

~2 T 


cosh 


Hiv^J (tanh A n ch) 


(A35) 


A steady-state form of the solution may be obtained by evaluating equation (A35) at V = 0. If this 
is done, the resulting equation for AT nmx is 


0((U) = Ar max 


&Lq o 

TT 2 ck r 


sin (mrw/2L) tanh ( nnch/2L ) 


(A36) 


Convergence of the Analytical Solution 

Since the transient and steady-state forms of the analytical solution (eqs. (A35) and (A36), 
respectively) are infinite series solutions, a convergence study was conducted to determine the number 
of terms necessary for an accurate solution. To simplify the convergence study, the steady-state 
analytical solution was used. However, the convergence criteria established as a result of the study 
were also applied to the transient analytical solution. The results of the convergence study are shown 
in figure 16. This figure shows the value of A T niax as a function of the number of terms used in the 
solution. As shown in the figure, the analytical solution oscillates about a converged solution. The 
solution is considered to be converged when the percent difference between successive maximum and 
minimum values is less than 1 percent. The percent difference is calculated with respect to the latest 
maximum. The temperature difference Ar[ nax is then calculated to be the average of the most recent 
maximum and minimum values. For the case shown in figure 16, convergence occurs at 320 terms. 
However, the number of terms required for convergence varies from case to case. 
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Effect of the Assumed Boundary Conditions 

To develop the analytical solution to the approximate problem, it is assumed that the tempera- 
tures along the edges £ = ±L are equal to T- m . Therefore, the value of 0 at these points is G°F. This 
boundary condition may produce erroneous results if it is placed too close to the square-heat-pulse 
region because it may act as a heat sink and artificially reduce the maximum temperature difference 
between the outer and inner surfaces. Therefore, the effect of varying the value of L on the maximum 
temperature differences predicted by the analytical solution was investigated using the steady-state 
analytical solution. The results of the investigation are presented in figure 17, wherein the converged 
value of AT max is shown as a function of L/w. (It is recognized that the temperatures shown in 
fig. 17 for the two composite materials exceed the maximum use temperatures for those materials; 
however, for the purpose of this discussion, the results are being retained.) As shown in the figure, 
the solution is unaffected by increases in L/w for L/w > 80. Thus, the value of L was set at 80tc in 
all analyses. 

Effect of Neglecting Radiation 

The effects of radiation heat transfer leaving the outer surface of the cowl leading edge are 
neglected in the analytical solution. To investigate the validity of this assumption, several cases are 
investigated with the steady-state finite-element analysis (with P-140/Incoloy 909) used while T[ n 
is varied from 0°F to 3000°F. The material P-140/Ineoloy 909 was selected for this investigation 
because it has the lowest values of transverse and longitudinal thermal conductivities; therefore, it 
can be expected to have the highest surface temperatures for any given conditions. The results of the 
investigation are shown in figure 18. As shown in the figure, the maximum temperature difference 
decreases from approximately 4950°F to approximately 1600° F as T Ui is varied from 0°F to 3000°F 
(a 7.1-percent decrease). The decreasing values of AT max indicate that radiation heat transfer 
leaving the outer surface of the leading edge does become more important as T lu is increased. Since 
the increased internal temperatures cause the maximum outer-surface temperatures to increase, the 
amount of heat being radiated away from the leading edge increases. However, because typical 
values of T\ n for a cowl- leading-edge structure should not exceed 2000° F, the effects of radiation 
heat transfer leaving the outer surface of the cowl leading edge art? neglected in the present study. 
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Table I. Material Properties 


(a) Constituent properties 


Material 

Btu/ft-hr-°F 

k T , 

Btu/ft-hr-°F 

p- 

lb/in* 

r i>< 

Btu/lb-°F 

Copper 

226 

226 

0.322 

0.10 

Incoioy 909 

9 

9 

.296 

.11 

P-140 carbon fiber 

450 

6 

.050 

.24 

HOPG graphite 

1156 

6 

.050 

.24 


(b) Composite properties (50-percent fiber- volume fraction) 


Material 

k L , 

Btu/ft-hr-°F 

k T . 

Btu/ft-hr-°F 

p- 

lb/in :! 

r-p. 

Btu/lb-°F 

P- 140/copper 

340 

12 

0.186 

0.17 

HOPG /copper 

690 

12 

.186 

.17 

P-140/Incoloy 909 

230 

7 

.173 

.17 

HOPG /Incoioy 909 

583 

7 

.173 

.17 


Table II. Ranges of Parameter Variations 


V , in/sec . . 
q G . Btu/ft 2 -scc 
ki, Btu/ft-hr-°F 
kq\ Btu/ft-hr-°F 
hi/ hr .... 

2uk in 

/i, in 


. . . . 0 to 20 

10 000 to 100 000 
. . . . 9 to 690 
. ... 7 to 226 
. . . 1.0 to 83.3 
. . 0.01 to 0.03 
. 0.005 to 0.030 
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Plume 
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reusable air-breathing hypersonic vehicle. 



l-leading-edge 
shock wave 



Figure 2. Schematic of Type IV shock-wave interference pattern (supersonic jet) impinging on a cylinder. 



Cowl-leading-edge bow 

shock wave Cowl-leading-edge bow 

\ shock wave -v 
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(a) Comparison of experimental and finite-element (b) Heating distribution assumed for transient case 
heating distributions. M « 8. (From ref. 2.) 



(c) Heating distribution assumed for steady-state case. 

Figure 4. Type IV shock- wave-interference heating modeled as uniform square heat pulse for transient and 
steady-state cases. 
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384 line heat source elements 
424 line radiation elements 
2688 degrees of freedom 


Figure 6. Detailed view of cylindrical portion of thermal-finite-elenient model of engine cowl leading edge 




Figure 8. Effect of square-heat-pulse sweep speed 
composites, h. = 0.02 in.; 2 w = 0.01 in.; q 0 = 50 


Finite-element Analytical 

analysis solution 

O P-1 40/copper 

□ HOPG/copper 

A Copper 


A Copper 

n 



A 


A 


A 


10 


15 


J 

20 


sp speed, V, in/sec 

maximum temperature differences for copper and its 
Btu/ft 2 -sec; T m = 0°F. 


Figure 9. Effect of square-heat-pulse sweep speed on maximum temperature differences for Incolov 909 and its 
composites, h = 0.02 in.; 2 w — 0.01 in.; qo = 50000 Btu/ft -sec, Tj n OF. 



F-irr-sec/Btu 



Figure 11. Effect of cowl- leading-edge thicl 
composites, w = 0.005 in. 



.010 .015 .020 .025 

Cowl-leading-edge thickness, h, in. 


is on maximum temperature difference for copper and its 





Cowl-leading-edge thickness, h, in. 


Figure 12. Effect of cowl-leading-edge thickness on maximum temperature difference for Incoloy 909 composites. 
w — 0.005 in. 
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Figure 14. Effect of square-heat-pulse half-width on maximum temperature difference for Incolov 909 
composites. 




Figure 16. Maximum temperature difference as fum 
solution. 
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Analytical results 
Converged solution 



240 320 400 

icluded in steady-state solution 


of number of terms used in steady-state analytical 




Solution unaffected by Uw increase 



Figure 17. Effect of position of lateral boundary condition used in analytical solution on maximum temperature 1 
difference, h = 0.02 in.; q 0 — 50 000 Btu/ft^-sec. 




Figure 18. Effect of radiation on maximum temperature difference of P-140/Incoloy 909 composite as T- m is in- 
creased. ki = 230 Btu/ft-hr-°F; k? — 7 Btu/ft-hr-°F; h = 0.02 in.; w — 0.005 in.; q Q = 50 000 Btu/ft 2 -sec 
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